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A unit evolution step of discrete-time quantum walks is determined by both a coin-flip operator 
and a position-shift operator. The behavior of quantum walkers after many steps delicately depends 
on the coin-flip operator and an initial condition of the walk. To get the behavior, a lot of long- 
time limit distributions for the quantum walks starting with a localized initial state have been 
derived. In the present paper, we compute limit distributions of a 2-state quantum walk with a 
non-localized initial state, not a localized initial state, and discuss how the walker depends on the 
coin-flip operator. The initial state induced from the Fourier series expansion, which is called the 
(q, /3) non-localized initial state in this paper, provides different limit density functions from the 
ones of the quantum walk with a localized initial state. 



I. INTRODUCTION 

Discrete-time 2-state quantum walks (QWs) are con- 
sidered as motions of quantum particles with the up-spin 
or down-spin state. The system of quantum walkers is 
described by probability amplitude, and its dynamics is 
determined by a unitary operator which is constructed 
of both a coin-flip operator and a position-shift operator. 
That is, the QW is a unitary process. One of the goals 
of the study is to know where the quantum walker is af- 
ter many steps. To find the answer, a lot of long-time 
limit theorems have been derived. Particularly, the limit 
distributions play an important role in letting us know 
the spatial distribution of the walkers after long time. 
While we theoretically investigate the QWs, experiments 
on the QW after some steps have been also realized in 
physical quantum systems recently. In preparation for 
the success of realization of the QWs after many steps, 
we require investigation of the limit theorems which can 
support justice of the experiments. 

The QWs have provided us with impressive limit dis- 
tributions. The limit distribution of the 2-state QWs on 
the line Z = {0, ±1, ±2, . . .} is depicted by a density func- 
tion with a compact support like the arcsine law [l|, y] . 
If the quantum walkers possess more than three states, 
they can spatially localize in probability distribution and 
the limit density functions are described by both a con- 
tinuous spectrum and the Dirac (5-function which means 
occasion of localization [^-Q ■ The limit density functions 
of the QW on the two-dimensional square lattice also 
have a compact support [T], l8|. Long-time asymptotic 
behavior of the QW on more than three-dimensional lat- 
tice is currently open. Outside lattices, the Grover walks, 
whose coin-shift operators are defined by symmetric uni- 
tary matrices, were analyzed on trees, joint half lines or 
a spidernet [9l-[ll|. These limit theorems are results for 
the QWs with a localized initial condition. The author. 



however, thinks why we focus on the QWs starting from 
a certain position always. To understand properties of 
the QWs, we need to clarify the relation between the 
QWs and more general initial conditions. In fact, when 
we study quantum systems, we often focus on initial con- 
ditions distributed over a domain. For example, in quan- 
tum search algorithms (e.g. Grover's search algorithm), 
we use quantum systems starting from a non-localized 
initial state. So, the aim of this paper is to analyze the 
QW on the line with a non-localized initial state, and 
we report contribution from the coin-flip operator to the 
limit distributions. When the walker on the line starts 
with a localized initial state, it is known that the coin- 
flip operator gives an effect to asymmetric shape of the 
limit distribution [l|, 0]. It would be interesting to ex- 
plore the dependence on the coin-flip operator in case of 
a non-localized initial state. 

Although the initial state given in this paper is special, 
our results would become the first step to catch the rela- 
tion between the limit distributions of the QWs and their 
non-localized initial states. By knowing how the QWs 
depend on the initial state, we can reach deep under- 
standing of their interesting characteristics. As a result, 
the QWs might become attractive in a field of applica- 
tion and spread into science more widely. There are a 
few results for the QWs distributed widely in an initial 
state dMll. Abal et. al. [ll[ll analyzed a QW start- 
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ing from two positions. Chandrashekar and Busch [15 1 
reported numerical results for a QW with an initial state 
ranging over an area. Valcarcel et al. [ij] treated the QW 
initialized by a Gaussian-like distribution in a continuum 
limit. 

The properties of the QWs have been clarified the- 
oretically, while there are some experimental results of 
the QWs, which were performed by trapped ions, atoms 
in optical lattices, photons, optical devices and Bose- 
Einstein condensates |16l - l24l |. Although a sizable num- 
ber of steps in random walks are easily realized, wc re- 
quire a hard effort to demonstrate many steps of quantum 
walkers. For example, Zahringer et al. [23| implemented 
QWs with up to 23 steps by using trapped ions. Since 



realization of the QWs is developing recently, we need to 
prepare rigorous results to contribute experiments on the 
QWs. 

The rest of this paper is organized as follows. In Sec. ITU 
we introduce the definition of a discrete-time 2-state QW 
on the line and the (a, /3) non- localized initial state which 
is a wider initial state. Section Hill is devoted to show our 
limit theorem. In Sec. lIV| we show contribution from the 
coin-flip operators to the limit theorem of the QW with 
the (a, /?) non-localized initial state in some examples. 
In the final section, we conclude our results and discuss 
a future problem. 



II. A DISCRETE-TIME QW ON THE LINE 
WITH A NON-LOCALIZED INITIAL STATE 

Total system of discrete-time 2-state QWs on the line 
is defined in a tensor space Hp ® He, where Up is 
called a position Hilbert space which is spanned by a 
basis {\x) : x £ Z} and He is called a coin Hilbert space 
which is spanned by a basis {|0) , |1)} with the vectors 
(0| = [1,0], (1| = [0,1]. Let \Mx)) e He be the prob- 
ability amplitudes of the walker at position x at time 
t e {0,1,2,...}. Physical meaning of the component 
(OjV'tla;)) (resp. {l\^ljt{x))) is the amplitude of the up- 
spin (resp. down-spin) state of the quantum particle at 
position X at time t. The state of the 2-state QWs on the 
line at time t is expressed by j^t) = J2x<£i, 1^) ® \^t{x)). 
Time evolution of our QWs is described by a unitary 
matrix, that is a coin-flip operator, 

U^^e =cos6'|O)(O| + (-l)«sin0|O)(l| 

+ sin6'|l)(0|-(-l)«cose'|l)(l| (1) 

with ^ e {0, 1} and 6 E [0, 27r). In this paper, we focus 
on the case oi 9 y^ Tr/2, tt, 3tt/2. We should note that 
Ui,^-g = (|0) (1| + |1) {0\)Uo,e anddetC/^,e = (-1)«. The 
amplitudes evolve along 



\^Pt+iix))=\0){0\U^^g\M=^ + l)) 



(2) 



Equation ([2]) means that the walker's position is shifted 
after its coin states are flipped. The probability that the 
quantum walker Xt can be found at position x at time t 
is defined by 



viXt = x) ^ {Mx)\Mx)) 



iMM^)) 



(3) 



The term | {0\iJt{x)) |^ (resp. | {l\Mx)) P) denotes the 
probability that the quantum particle with the up-spin 
(resp. down-spin) state can be found at position x at 
time t. By giving the initial states \^poix)), we get the 
probability distribution P{Xt = x) for any time t. In this 
paper, we focus on a special non-localized initial state as 
follows. 



Let wi,W2 '■ IR — > R be the functions that satisfy 
Wj{k + 2tt) = Wj{k) {j = 1,2) and Wj S L^[— tt, tt]. As- 
suming J^^ {\wi{k)\'^ + \w2{k)\'^] dk > 0, we build the 
following initial condition of the 2-statc QWs, 



1^-0(2;)) = 



y^W{wi,W2) 



{d^{x)+td2{x)}\cj,), (4) 



where 



W{WI,W2)= {\wi{k)\^ + \w2{k)\^}dk, (5) 

J —7T 

1 



dj{x) 



V^TT J-TT 



w,{k)e'''^dk (j = l,2), (6) 



and 

l«|2 



\^) 



ajO) -I- ;3|1) with a,^ e C and jaj^ -|- 
1. We should note that Y.xez'^i^'^ = x) = 
X^aiez ('/'o(^)IV'o(a;)) = 1 by Parseval's theorem. Since 
Eq. (|4]) provides a non-localized initial state, we call it 
the (a, /3) non-localized initial state in this paper. 



HI. LIMIT DISTRIBUTION 

Limit distributions depict approximate motion of the 
QW after many steps, just like the relation between ran- 
dom walks and central limit theorems. To get the be- 
havior of the QW with the (a, /3) non-localized initial 
state as t — >■ 00, we compute a limit distribution by 
using the Fourier analysis which is one of the standard 
methods to derive limit theorems of the QWs 0-H, l25l - 
Wf\ . From Eq. ^, the time evolution of the Fourier 
transform \i>t{k)) = Y.x&^'''''" \^t{x)) {k e [-7r,n)) 
becomes \i>t+i{k)) = U^,e{k)\^tik)) = U^MkY\^o{k)) 
, where U^^e{k) = {e^^ |0) (0| + e-^^ |1) {l\)U^fi. The 
Carleson-Hunt theorem changes the initial state of the 
Fourier transform. 



*o(fc))=E (^/ ^liky'^'dlAe 



X&L 



.,±/>.(iK-.i)e-L/^^J;^ 



27r 



W[W\,W2) 



{u;i(A:)-fW2(fc) }!</'), 



(7) 



almost everywhere on [— 7r,7r] [2^, [2^. Equation ([7]) is 
essentially the Fourier series expansions of the functions 
wi{k),W2{k). After using the method treated in Grim- 



mett et al. [27|, for r = 0, 1, 2, . . ., we get 



lim E 



t 



rlh 



(i;(fc)|J4l^o(fc + W2)) 



(«(-fc)|J^|*o(-fc + evr/2)) 



+ (-ir<j («(7r-fc)|J4|vI/o(fc + W2)> 

_ . 2 " 



(«(7r + A:)|J^|*o(-fc + W2)) 



where 
/i(fc) 

\v{k)) = 



coos fc 



V 1 — (? sin fc 



(8) 



(9) 



e'^^s 



v/]V(fc) 






N(k) =1 + s^ + (? cos 2fc + 2c cos fc Vl - c^ sin^ fc, 

(11) 
Je=|0)(0| + (-l)«|l)(l|, (12) 

and c = COS0, s = sinfl. By putting /i(fc) = x, we have 



hm E 

i— j-oo 






> 27r(l - x2)Vc2 - a;2 

X ry(a;)/(_|c|,|c|)(2;)c^a;, (13) 
where /^ (a;) = 1 if a; e A, /^ (a;) = if a; ^ A and 

77(a;) = |(«(«;(a:))| J^ |#o(a^(.t) + W2))f 

(«(^(a;))|J4l^o(«(a:)-7r + evr/2))|' 

(«(-K(x))|J^|^o("«;(a:)+evr/2))|' 

+ I (t^(-Aj(x))| J^ |^o(^ - n{x) + eV2)) f , (14) 
|s|a; 



At (a;) — arccos 



z^f\ 



€[0,7r]. 



(15) 



In Eq. (|13p , if we give a special initial state to the Fourier 
transform |v['o(fc)), we can calculate a limit distribution. 
Let F : R — > C be a function that satisfies F(k + 27r) — 
F{k), J^^\F{k)\'^dk = 2n and 3?(F(fc)), 9(F(fc)) S 
C°°[— TT, tt] almost everywhere, where ^{z) (resp. 3(2)) 
denotes the real (resp. imaginary) part of z e C and 
R {resp. C) means the set of real (resp. complex) num- 
bers. If we set l^o(fc)) = P{k) \(t>}, we obtain 



hm E 



t 



x'' lfiix;^,a,l3)r]i{x;^) 
+ f2ix; ^, a, P)r]2ix; ^)j I(^^\c\,\c\)ix) dx, (16) 



where 

/i(x;^,a,/3) 



7r(l — a;2)vc? 

|2 



1- 



/2(x;e,a,/3)=-(-l)^ 



, ^ ^ ^^^ 2gJft(a/3) 

c 
s5(a;5) 



(17) 
(18) 



(19) 



|c|7r(l-a;2)' 

ryi(x;e) = J{|F(«;(a;)+eV2)P 

+ |F(-«(a;)+eV2)P 
+ |F(K(a;)-7r + e7r/2)|^ 
+ |F(^-/c(a;)+eV2)|2} 

^2(x;0 = ^{in^(:^)+W2)P 

-\F{-Kix)+^7r/2)\^ 
+ \F{k{x)~tt + £,7t/2)\^ 

-\F{TT~K{x)+i7:/2)\'] 

As you can see from Eq. ([7]), the function F{k) corre- 
sponds to 



(20) 



w(wlw2) {^i'i(fc) + iw2{k)}. In the next sec- 
tion, we will show contribution from the coin-flip opera- 
tors U(^^e to the limit theorem Eq. ^TE\\ when the functions 
wi{k),W2{k) are supplied. 



IV. CONTRIBUTION FROM COIN-FLIP 
OPERATORS 

In this section, to see contribution from the coin- 
flip operators f/{,6/, we concretely give the functions 
wi{k), W2{k) on a range with the length 27r and compute 
limit density functions of the QW as i — ^ 00. Outside 
the given domains of wi(fc), W2{k), the functions shall be 
27r-periodically expanded on R. 



A. Case 1 

We take wi{k) = cosafe(— tt < k < 7r),W2{k) = 
sinafc(— TT < fc < 7r) with a ^ 1. Then the (a,/3) non- 
localized initial state becomes 



\Mx)) = (-1)^ 



Sm flTT 

(x -I- a)7r 



(a; = 0,±l,±2. 



(21) 

Before we show the limit distribution for this initial state, 
we remark the following fact from Eq. ([T6| . If the func- 
tion F(fc) given in SecinDsatisfies |F(fc-7r)| = \F{-k)\ ^ 
|i^(fc)|, then we have 



lim E 



t 



x" h(x\i,a,Pi) 



X \F{k{x) + ^7r/2)|2/(_|,|_|,|)(a;) dx. (22) 



By using this result, we compute the hmit distribution 
of the QW with the initial state Eq. (pij) . 



lim E 

t— >oo 


IfXtV' 

[\t J \ 



x^'fiix; £,, a, /3)/(_|c|,|c|) {x)dx, 



■■I{-\c\^\c\){x) 



7r(l — x'^)\/c'^ — x^ 

'I - l/3| +(-1)'" —^ >x 



dx. (23) 



Equation ([23]) is the same limit theorem as that of the 
walker starting from the origin with |V'o(0)) = \4>) [H, 0I- 



-1 -|c| 



1 -1 -|cl 



|c| 1 



(a) C = 



(b) e = 1 



FIG. 1. (Case 1) Comparison between the limit density func- 
tion (red line) and the probability distribution f{Xt/t — x) 
at time t = 5000 (blue line) with a = 1/^2, /3 = i/V2. 



B. Case 2 

Given the functions wi{k) — log(4cos^fc) {tt/2 < k < 
57r/2), W2{k) = 0, the initial state is expressed as 



\M^)) = 



01 



'(-I)- 



. 2V3j£| 



(a; = 0,±l,±3,...), 
(x = ±2,±4,...). 



In a similar fashion as Case 1, we can extract 



(24) 



hm E 


\(XtV] 

[\t J \ 



x''fi{x;^,a,l3) 

D 

'X gi{x;Ol{-\c\,\c\){x)dx, (25) 



where 



3 f is^x^ 
9i{x; 0) =— < log ^- 5- 

, , 3 r, 4(c2~2;2) 
31 (a:; 1 =— <^ log ^- 5- 

TT^ [^ C^(l — X^) 



(26) 

(27) 



We remark that the limit density function with ^ = has 
three singular points x = 0, ±|c|. 



C. Case 3 

The (a, /?) non-localized initial state constructed 
of wi{k) = ilog(cot|A:/2-7r/4|)(-7r/2 < fc < 



\r\ 1 -1 -IcI 



(a) C = 



(b) C = 1 



FIG. 2. (Case 2) Comparison between the limit density func- 
tion (red line) and the probability distribution f{Xt/t = x) 
at time t = 5000 (blue line) with a = l/\/2, /3 = i/\/2. 



37r/2), W2{k) = is written as 



\M^)) = 



o|</.) 



■ 2 



(x = 0,±2,...), 
(x = ±l,±3,...) 



(28) 



The functions wi{k),W2{k) don't satisfy the condition 
corresponding to |F(fc— 7r)| = \F{—k)\ — \F{k)\. We have 
to pass through the different computation from Eq. (|22l) . 
So, we add a result derived immediately from Eq. p^ as 
follows. If the function F{k) satisfies \F{-k)\'^ + |F(7r - 
fc)|2 ^ |F(fc)|2 + \F{k - 7r)|2, then we have 



lim E 



t 



x''fi{x;^,a,(3) 
'xV3{x;0l{-\c\M){x)dx, (29) 



where 



+ \F{K{x)-n + £,TT/2)\^y (30) 



From Eq. (I29p . we obtain the quantum walker with 
Eq. (US]) as t ^ 00, 



lim E 



where 
g2{x;0) 



t 



x''fiix;^,a,f3) 
'x- 92{x;0l{-\c\,\c\){x)dx. (31) 



log 



:|vr 



-|- \/c^ — x"^ 



+ <log 



g2{x;l) = -ir 



log 



lc|v^r^^^-%/?" 



(32) 



log 



c|c|Vl — x^ + \cs\x 



c\/ (?■ — x^ 



c\c\\/\ — x"^ — \cs\x 
c\/ (? — x"^ 



2-1 



(33) 



Note that the limit density function with ^ = isn't 
continuous at a; = 0. On the other hand, since 
lim2,_>±|c|/i(a;;0, a,/3)52(a;;0)/(_|c|,|c|)(a;) = 0, it's con- 
tinuous at X = ±|c|. 




(a)C 



(b)C^ 



L 



(a) C = 



(b)C 



FIG. 3. (Case 3) Comparison between the limit density func- 
tion (red line) and the probability distribution f(Xt/t — x) 
at time t = 5000 (blue line) with a = l/\/2, P = i/V2. 



FIG. 4. (Case 4) Comparison between the limit density func- 
tion (red line) and the probability distribution f{Xt/t = x) 
at time t = 5000 (blue line) with a = l/V^, P = i/V2. 



D. Case 4 



E. Case 5 



Let us suppose 

{0 (-TT < A: < -7r/2), 

sinfc+1 (-7r/2 < fc < 7r/2), (34) 

2 (7r/2 < A: < tt), 
W2{k) =0. (35) 

These functions yield the following (a,/3) non- localized 
initial state, 



\Mx)) = < 






Equation p^ has the representation 



{x = 0), 

(a; = ±l), 

(x = ±3, ±5, . . .), 

(a; = ±2, ±6,...), 
(a; = ±4, ±8,...). 

(36) 



lim I 



where 



t 



x''{fi{x;ta,P)g3{x;C) 



+ h^x; C, a, /3)g4(a;; C)}/(-|c|,|c|)(a;) dx, (37) 



(38) 
(39) 



53(.;0) = -|3 + -,^^— ^ 

2 f s^a;^ 



54(2;; 0) 



54(2;; 1) 



(1 



8c 



\\% 



\J c^ — x^ 



_Jsja 



(-|c|<x<0), 



7 V|c| 



=\/T^ 



(0<x<|c|), 
(-|c| <x<0). 



(40) 



(41) 



Notice that the limit density function with ^ = 1 isn't 
continuous at a; = 0. 



As the final example, we focus on an initial state over 
a finite domain. For a non- negative integer n, we give 

u;i (fc) = cos2"+2 fc -f sin2"+2 fc, (42) 

W2(fc) =(-l)"(cosfcsin2"+ifc--sinfccos2"+ifc). (43) 

Then setting Xn^i = (-1)"+^ {(-lp(2j + 1) - l} (j = 
0, 1, 2, . . . ,n), we have the (a,/3) non-localized initial 
state, 



IV-oCx)) 



1 /v^r 

2^"Y 2r(2 

0|(/.) 


2n+2) /2n 
n+l) Vn 


(Vj, 


a; = 
x + 


Xr 
Xr 


(44) 



°° oZ-1. 



' ds is the Gamma function. A 



where T{z) — Jn s 

limit distribution follows from Eq. (|22|) : 

X 55(x)/(_|c|,|c|)(a;)c?a;, (45) 



lim E 



^ 

f 



where 

55(2;) = 



0Fr(2n + 2) s4n+2^4„+2 ^ (^2 „ 2;2)2«+l 



2r(2n 



l^4n+2(l _ j.2\2n+l 



(46) 



If n = 0, then the limit distribution is equivalent to that 
of the QW starting from the origin [l|, y] . The effect of 
the parameters ^ S {0, 1} is expressed in a linear term of 

/i(x;^,a,/3). 



V. SUMMARY 

In this section, we refer to our results and a future 
problem. Limit theorems of QWs have been mainly in- 
vestigated for localized initial conditions. Theoretical 
analysis for the QWs with a non-localized initial state 
is generally complex. In this paper, to deepen our under- 
standing of relationship between a 2-state QW on the 
line and its initial conditions, we produced the (a, /3) 
non-localized initial state induced from the Fourier se- 
ries expansion, and derived a long-time limit theorem for 



I 0.02 

1 



^^ 



-1 ~|c] 

X 

(a) ¥{Xt/t = x) 



(h) Mx;^,a,P)g,{x) 



FIG. 5. (Case 5) Comparison between the probability distri- 
bution ¥{Xt/t = a;) at time t = 1000 (Fig. (a)) and the hmit 
density function (Fig. (b)) with n — 50 and a = l/-\/2, (3 = 
i/V2. 



the walk. That is, we found one of the QWs with a non- 
locaUzed initial state whose approximate behavior after 
many steps can be forecasted. By using our limit the- 
orem, we have specifically computed limit density func- 
tions for some QWs with the (a, /?) non-localized initial 
state. As a result, Eqs. ^, dSJ), dST]) and (gS]) are 
different from the limit distribution for the 2-state QW 
starting with a localized initial state. The limit distribu- 



tions, moreover, taught us the effects of the coin-sift oper- 
ators U^ g = cos6l |0) (0| + (-!)« sine* |0) (l|+sin6' |1) (0|- 
(-l)«cos6'|l>(l| (ee{0,l}) (See Figs. IE] and H. 

As a development of our investigation, it is interesting 
to analyze long-time behavior of multi-state QWs with 
a non-localized initial state. Since the multi-state QWs 
starting from the origin can localize in probability dis- 
tributions after long time [3|-|5|, it is a conceivable topic 
to clarify the difference between our 2-state QW and the 
multi-state QWs with a non-localized initial state. Also, 
our result claims that we can create various probability 
distributions with a compact support from the 2-state 
QWs by controlling the functions wi{k),W2{k). If the 
QW after many steps is realized in a quantum system, 
we'll physically gain our desired probability distributions. 
At the same time, that means we can obtain our desired 
random variables in quantum computing. The author 
hopes that our results will contribute to the quantum 
computation theory. 

The author acknowledges support from the Meiji Uni- 
versity Global COE Program "Formation and Develop- 
ment of Mathematical Sciences Based on Modeling and 
Analysis" . 
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